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We give a variational formulation of classical statistical mechanics where the one-body density and 
the local entropy distribution constitute the trial fields. Using Levy's constrained search method 
it is shown that the grand potential is a functional of both distributions, that it is minimal in 
equilibrium, and that the minimizing fields are those at equilibrium. The functional splits into a sum 
of entropic, external energetic and internal energetic contributions. Several common approximate 
Helmholtz free energy density functionals, such as the Rosenfeld fundamental measure theory for 
hard sphere mixtures, are transformed to internal energy functionals. The variational derivatives of 
the internal energy functional are used to generalize dynamical density functional theory to include 
the dynamics of the microscopic entropy distribution, as is relevant for studying heat transport and 
thermal diffusion. 
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I. INTRODUCTION 

The study of classical many-body systems in equilib- 
rium is often based on the grand potential D,q as a func- 
tion of its natural variables, which for a one-component 
system are the chemical potential p, the temperature T 
and the volume V. The mean number of particles, Nq, 
is then obtained as a partial derivative dflo/dfi = —No, 
while keeping T and V fixed, and the mean particle den- 
sity is simply Nq/V. When the Hamiltonian contains a 
contribution due to an external potential, v(r), where r is 
the position coordinate, then the density is in general no 
longer uniform, but becomes position-dependent, hence 
Po{r), where po(r)dr gives the mean number of particles 
in a volume element dr. Here the difference p — v(r) 
and po( r ) play the role of conjugate thermodynamic 
fields, and for convenience one often defines (formally) 
a position-dependent chemical potential p(r) = p — v(r). 
The one-body density distribution can then be obtained 
as a functional derivative 8ilo/8p(r) = ~pa(r). 

Density functional theory (DFT) |l[ amounts to gen- 
eralizing this concept to a functional dependence of the 
grand potential on the one-body density distribution, 
i.e., going from Oq(//, T, V) to a functional Q([p],fi, T, V), 
where pM is a trial field. The variational principle of 
DFT [lj,l2| states that for given thermodynamic state (i.e. 
fixed values of T, p, and V) the density distribution that 
minimizes f2 is the physically realized equilibrium density 
po( r )- The non-trivial (additive) contribution to fi is the 
Helmholtz free energy functional F([p],T, V), which is in- 
dependent of p and generalizes the equilibrium Helmholtz 
free energy Fo(No,T,V) to a functional dependence on 
the trial density p(r). Inserting the equilibrium den- 
sity into the functional yields the equilibrium free energy, 
i.e., F (N , T, V) = F([p Q ],T,V), where N = J dr Po (r). 
There is a significant body of literature on application 
of this framework to a wide variety of interesting many- 
body phenomena in liquids and solids [3|-|5|- Both con- 
ceptually, and in practical DFT applications, the tem- 
perature enters as a mere parameter, often in the form 



of "thermal energy" fc^T, where ks is the Boltzmann 
constant. Clearly this situation is very different from 
the sophisticated treatment that the chemical potential 
received via introduction of p(r) and its conjugate field 
p(r). One might justify this by the fact that p(r) ^ const 
leads to a well-defined equilibrium when v(r) acts on the 
system, whereas considering inhomogeneous temperature 
distributions reeks of non-equilibrium. 

In thermodynamics, one can proceed and Legendre 
transform to the internal energy Eq(Nq, S,V), where 
the entropy S is an extensive state variable, conjugate 
to temperature T. The latter is recovered from T = 
dE /dS at Nq,V = const. Were one to generalize to an 
internal energy functional, one needed to introduce (and 
define) an entropy distribution s(r) that would "localize" 
(i.e. make dependent on position) the bulk entropy per 
unit volume, S/V. This programme possesses several 
requirements in order to be rigorous, i) A microscopic 
definition of the entropy distribution s(r) needs to be 
given, ii) The grand potential functional needs to depend 
on both the microscopic density and the microscopic en- 
tropy, i.e., Q([p,s],p,T,V). Its non-trivial contribution 
should be an internal energy functional of both micro- 
scopic distributions, i.e., E([p,s],V). iii) The general- 
ized grand potential functional should be minimal at the 
equilibrium values po( r ) an d s (r). Note that T — const 
will be associated in general with a non-trivial spatial 
dependence so( r ) 7^ const. This forms a generalization 
of the simple parametric dependence on temperature in 
conventional DFT to a proper Euler-Lagrange equation. 

In the following such a framework is established. We 
formulate the variational principle in Sec. |Hj Standard 
DFT approximation are converted to the internal energy 
representation in Sec. IIIII This includes internal energy 
functionals for the ideal gas, hard spheres in the funda- 
mental measures approximation d i-i], the quadratic 
mean-field functional Q, etc. Based on the continuity 
equations for particle density and internal energy den- 
sity, and inspired by the framework of linear irreversible 
thermodynamics, in Sec. IIVI we formulate a theory for 
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diffusive dynamics that corresponds to dynamical DFT 
(DDFT) 0, US [n|, but includes the dynamics of the 
entropy current. Conclusions are given in Sec. |V] 



II. VARIATIONAL PRINCIPLE 

We consider a classical system with N particles and 
Hamiltonian Hpj. The equilibrium many-body probabil- 
ity distribution in the grand ensemble is given by 



fo = 3 1 exp - 



H N - pN 
k B T 



(1) 



Here the normalization constant is the grand partition 
sum 



3 = Tr exp 



H N - (j,N 
knT 



with the (classical) trace being defined as 
1 



Tr 



^ WjvT f dri --- dTN fdpi...dp N , 

N=0 ' J J 



(2) 



(3) 



where h is Planck's constant, is the position coordi- 
nate and Pi is the momentum of particle i = 1,...,N. 
Mermin's form Q for the grand potential as a functional 
of a trial many-body distribution / is 



Sl[f]=Trf(H N - f iN + k B Tlnf). 



(4) 



Here / is an arbitrary many-body distribution that is 
normalized, i.e., that satisfies 

Tr/ = 1. (5) 

Inserting the equilibrium distribution ([TJ into (j4} yields 



fi[/o] = Tl ' fo { h n - »N + k B T\n f ) 



(6) 



Tr/o 



H N -»N-k B T(ln3+^^t 



knT 



= — fcsTlna 



(7) 
(8) 
(9) 



where f^o is the equilibrium grand potential. From the 
Gibbs inequality it is straightforward to show [l|, [l^] that 
for any / ^ / the inequality fi[/] > fi[/o] holds and 
hence 

fi = minTr f(H N - pN + k B T\nf), (10) 

We use the conventional definition of the density op- 
erator 



N 



(11) 



where £(-) is the (three-dimensional) Dirac distribution, 
and express the one-body density distribution in equilib- 
rium as the average 



P o(r) =Trp(r)/ . 



(12) 



We also define a position-dependent entropy density (per 
unit volume) as 



So (r) = -fc B Tr^/ ln/ . 



(13) 



Note that the integral -T J drs a (r) = Trfc B T/ ln/o 
equals the entropic contribution to the grand potential, 
cf. the last term in Eq. 

We use Levy's constrained search method [H, [Hl|, as 
proved useful for classical systems [TEj], and express (fT0|) 
as a two-stage minimization 

O = min min Tr f(H N - pN + k B T\n /), (14) 

p,s f^-p,s 

where the inner minimization is performed for all trial / 
under the constraint that these generate the given density 
distribution p(r) and the given local entropy distribution 
s(r) via 



p(r) =Tr/5(r)/, 

s(r) = -fc s Tr^/ln/. 



(15) 
(16) 



The relationships (|T5|) and (fTo| are indicated as / — >• p, s 
in the notation of (p~4|) . 

In the following we restrict ourselves to Hamiltonians 
that consist of kinetic energy and internal and external 
contributions to the potential energy, i.e., that are of the 
form 



H 



N 



N 2 N 

E£ + C/ (ri>---^)+E^ 



2=1 



(17) 



where pj = p; ■ p,, m is the particle mass, U is the 
interparticle interaction potential, and v(r) is an external 
potential that acts on the system. Hence (fT4")l is more 
explicitly 



N o 

Pi 



r» = min min Tr/ Vf + (/(r 1 ,..,r I ) 
p.s f-tp,s \ * — ' Ira 
\i=\ 

N \ 

+ J2 v (n) ~ »N + k B T In/ J. (18) 

i=l / 

In the expression above several contributions can be writ- 
ten as space integrals over averaged one-body quantities. 
First, the terms due to the external and the chemical 
potential are 

Tr / (j>( r *) - ^ = J drp(r)(v(r) - M ), (19) 
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because / — > p via (fl5|) . Furthermore, the last term in 
(Ha is 



Tr fk B T In / = -T / drs(r), 



(20) 



because (|T6|) implies that / — > s. Hence the terms (|19|) 
and (|2"0"|) are constants with respect to the inner mini- 
mization in (1141) . Hence we can separate them out and 
arrive at 

O = mm ^E\p, s) + J dr [p(r)(u(r) - M ) - Ts(r)]| , 

(21) 

where we have defined the internal energy as a functional 
of the density and entropy distributions as 



E[p. s] = min 



N 



(22) 



where, once more, the minimization ("search" [13j) is 
constrained to all trial / that generate the given p(r) 
and s(r) via (fl~5|) and (fl"6|) . respectively. Here and in the 
following we suppress the dependence on volume V in the 
notation. 

The grand potential functional is the object inside of 
the minimization in (|21[) , defined as 

Q([p, sip, T) = E[p, s]-T J drs(r) + / drp(r)(v(r) - p). 

(23) 



Eq. (|2"Tj) then becomes 



£1 = minO([p, s],p,T), 



(24) 



which implies that the following functional derivatives 
vanish at equilibrium 



sn([ P ,s],p,T) 



Mr) 



= and 



6n([ P , s ],p,T) 



6s(r) 



= 0. 



(25) 



The density and entropy distribution that satisfy (|25[) 
are indeed po( r ) and So(r), as can be seen from their 
definitions, (fT2j) and ()13|) . and the fact that /o minimizes 
f2[/]. This implies that 

n (fi,T) = n([p , So },fx,T) (26) 

= E[p Q ,s ]-T J drs (r)+ / drp (r)(w(r) - p), 

(27) 

and that the internal energy in equilibrium is 

E (N ,S )=E[p ,s ], (28) 
where So = J drs {r). 



Using the definition (l23l) the Euler-Lagrange equations 
(|23|) can be cast in the form 



6E\p,; 



6p(r) 
5E[p,s] 



5s(r) 



Po,so 



Po,so 



p - v(r), 



T. 



(29) 
(30) 



For completeness, the Helmholtz free energy functional, 
on which DFT is conventionally built, is obtained as 

F([p],T) = mm \E[p, s]-T J dvs(v)^ (31) 



= E[p, s p ] - T J drs p (i 



(32) 



where s p (r) denotes the entropy distribution at the min- 
imum in (1311). which hence satisfies 



6E\p,s] 



5s{r) 



T. 



(33) 



where p(r) is the (trial) density distribution on the left 
hand side of ([31~|) . 

Eqs. ([2TR) and (|3U1) constitute a closed system of equa- 
tions for the determination of /9o(r) and so(r) for given 
thermodynamic statepoint /i, T and given external po- 
tential v(r). In practical applications one is required to 
use an approximation for E[p,s\. Hence it is interesting 
to formulate common free energy DFT approximations in 
the internal energy picture, as we do in the next section. 



III. EXAMPLES FOR INTERNAL ENERGY 
FUNCTION ALS 

We start with the ideal gas, where U(yi, . . . , rjv) = 0. 
The Helmholtz free energy functional can be derived from 
the fact that the absence of interactions decouples all vol- 
ume elements of the system P, QjJ ■ In each volume ele- 
ment the (bulk) ideal gas properties holds. Hence the free 
energy functional is an integral over a local free energy 
density, 

F id ([p],T) = k B T J dvp{v) [ln(p(r)A 3 ) - l] , (34) 

where the thermal de Broglie wavelength depends on T 
and is given by 



A = 



/ 2nh 2 

rnksT ' 



(35) 



and h = h/(2-K). Corresponding reasoning leads to the 
internal energy functional either by starting directly from 
the expression for the bulk internal energy of the ideal 
gas, or by Legendre transforming the integrand in (|34p . 
One arrives at the result 



Ei d [p,s] = 4S- / drp{rf/ 3 exp 

e ' 777. J 



Mr) ' 
3k B p(r), 



(36) 
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where e is the exponential constant. This result is the 
same as that obtained from Legendre transforming each 
volume element. The functional (|36|) is local and non- 
linear. Note that H 2 /m carries units of energy x length 2 , 
as is consistent with the integrand that has units of 
length -5 . Clearly, the comparison of (|36|) to (|34| points 
to the quite striking density power of 5/3 in ()36[) . and 
the fact that the entropy density per unit volume, s(r), 
appears in effect as an entropy density per particle, 
s(r)/p(r). 

Evaluating the derivatives in the Euler-Lagrange equa- 
tions (j2l)j) and (|3H)l and rearranging yields 



PoO) = A d exp - 



M — v ( r ) 
knT 



(37) 



«o (r) = ~k BPo (r) [ln(p (r)A 3 ) - 5/2] , (38) 



and insertion into (|36|) gives the internal energy of the 
ideal gas, solely due to kinetic contributions, Ei d [po, s o] = 
3k B T J drp(r)/2, a result which is certainly as expected. 

For interacting systems the total Hclmholtz free energy 
is usually split into an ideal and an excess (over ideal) 
contribution as 



F([p],T) = F id ([p},T) + F oxc ([p],T), 



(39) 



where F id ([p],T) is given by ([34]) and F exc ([p],T) de- 
scribes the effects of interparticle interactions. For hard 
spheres, most approximate functionals can be written in 
the form 



F exc ([p],T) = k B T / dr$(r), 



(40) 



where $(r) is a scaled excess free energy density per 
unit volume, which is independent of T . Temperature 
enters only via the global scaling factor k B T. For non- 
local functionals $(r) is a functional of p(r), typically via 
convolution. When such additional convolution integrals 
are present in the functional the choice which integral 
features as the "outer" integral in (|40p is not necessar- 
ily unique; see appendix [X] for a discussion of a suitable 
choice in fundamental-measure theory 0, 0] ■ Consider 
the following form of the internal energy functional 



#HS [P, S] = 



37rfr 2 
e 5 / 3 77 



drp(r) 5/3 exp 



K r ) - SHs([p],r) 



3fc s p(r)/2 



(41) 



where shsQp], 1 ") = — fe.B$(r) is the hard sphere contri- 
bution to the entropy. Eq. (j4Tj) is equivalent to (|40|) as 
can be seen from evaluating the Euler-Lagrange equa- 



tions ([291) and (|30|). which yield 
s (r) = -k BPo (r) [ln(p (r)A 3 ) - 5/2] + s HS ([po], r)], 



(42) 
(43) 



where c^g is equivalent to the one-body direct correlation 
function for hard spheres and obtained here as 



<&kM.r) = ^J^) j dr'snsdpW). (44) 

The common random phase approximation (RPA) [l2j 
consists of splitting a given interparticle pair interaction 
potential <p(r), where r is the particle-particle distance 
and J7(ri,...,rjv) = X)i<i<Kl r i _ r jl)> into a short- 
ranged repulsive, say hard sphere part 4>ns(r) ancl a long- 
ranged and slowly varying contribution <fi*(r), so that 
4>(r) = 0Hs( r ) + 4>*( r )- The corresponding internal en- 
ergy functional is 

E RPA [p,s}=E HS [p,s} + ^ J dvj dr'p(r)p(r')M\r~r'\), 

(45) 

where the effects of (j)tis(r) are described by the hard 
sphere functional (j4Tj) . Eq. (j45|) leads to the same "en- 
tropic" Euler-Lagrange equation (14"2"T) as for hard spheres, 
because 6E R pA/5s(r) = 5Ens/5s(r), and generates an 
additional contribution — J dr'p(r')^*(|r — r'\)/(k B T) in- 
side of the exponential in the "density" Euler-Lagrange 
equation (|4"3")) . 

For any system where the bulk internal energy 
Eo(No, S,V) is known, division by volume yields an in- 
ternal energy density e (p, s) = E (N /V, S/V)/V, from 
which in a local density approximation (LDA) the func- 
tional 



£lda[p, s] = J dre (p(r), s(r)) 



(46) 



follows. This is expected to be a good approximations 
when the smallest length scale over which p(r) and s(r) 
vary is much larger than all correlation lengths in the 
system. 

A further "generic" approximation, analogous to the 
Ramakrishnan-Youssouf (RY) [l6j quadratic approxima- 
tion, is to truncate the functional Taylor expansion 
around a homogeneous state with p(r) = pb = const and 
sir) = s b = const at second order in density, 



E RY [p,s] = E b (p b ,s b ) + ^ J dvj dr' [Ap(r)Ap{r')b pp (\r - r'|) 

+2Ap(r)As(r')M|r - r'|) + As(r)As(r')M|r - r'|)] 



(47) 
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where Ap(r) = p(r) — p& and As(r) = s(r) — Sb are the 
deviations from the respective bulk values, the subscript b 
indicates bulk quantities and the b pp {r), b ps (r) and b SB {r) 
are the second functional derivatives of E[p, s] evaluated 
in the homogeneous bulk, 



M|r-r'|) = 



8 2 E[p, 



5a(r)Sb(r') 



In general, 



b PP (r,r') = 
b ps (r,r') = 
b ss (r,r') = 



6 2 E{p,; 



8p{v)5 P (v>) 
S 2 E[ PlS ] 



5p(r)Ss(r') 
S 2 E[p,s] 



6s(r)6s(r') 



a,b = p,s. (48) 

(49) 
(50) 

Po.so 

(51) 



Note that the first order terms in (|4"T)) vanish, as one ex- 
pands around equilibrium and hence the Euler-Lagrange 
equations (|2T)1) and (13U)) hold (in the case v(r) — 0). 

Note that (|49|) - (f5"Tj) are analogous to the usual two- 
body direct correlation function obtained from the excess 
Helmholtz free energy functional as 



c 2 (r,r') = -(k B T) 



_i S 2 F cxc ([p],T) 



5p(r)5p(r') 



(52) 



IV. DIFFUSIVE DYNAMICS 

Using the equilibrium framework developed in Sec. HH 
we find it interesting to use it in a dynamical context, 
similar in spirit to dynamical density functional the- 
ory (DDFT) which rests on the equilibrium Helmholtz 
free energy functional. Much current research activity 
is aimed at applying and developing DDFT, which pro- 
vides a dynamical equation for the time evolution of the 
density profile p(r,i), where t is time. In order to derive 
such an equation the continuity equation for the density 
profile, which is exact, is supplemented by approxima- 
tions for the "thermodynamic driving force" that acts on 
the density. As compared to a diffusion equation gradi- 
ents in chemical potential are replaced by gradients in the 
functional (density) derivative of the Helmholtz free en- 
ergy density. Starting from a more microscopic point of 
view, DDFT can also be derived from the Smoluchowski 
equation [loj . 

Here, we spell out a similar framework for the joint 
time evolution of p(r, t) and the time and position- 
dependent entropy distribution s(r, t). We keep the dis- 
cussion at a phenomenological level and make no at- 
tempts at a derivation from first principles, albeit pay- 
ing attention that fundamental symmetry relations, i.e., 
the Onsager reciprocal relations, are satisfied. Hence the 
strategy consists of taking the appropriate dynamic equa- 
tions from linear irreversible thermodynamics 



ic equa- 

dam 



and replacing the temperature and density fields in the 
continuum description by the microscopic (functional) 
derivatives of the internal energy functional. 

We impose two continuity equations, one for the den- 
sity p(r, t) and one for the internal energy density e(r, t), 



p(T,t) = -V'3 p (T,t), 

e(r,f) = -V-J e (r,t), 



(53) 
(54) 



where the dot denotes a partial time derivative, i.e. 
p = dp/dt and e = de/dt. Solving the Gibbs-Duhcm 
relation de = Tds + pdp for the differential entropy per 
unit volume, ds, gives 



ds 



^de - j-dp, 



(55) 



from which the prefactors of the differentials on the right 
hand side are identified as the "driving forces" for the 
internal energy current J e and for the particle density 



current J p . Hence 



J p = DpV- 



-fi 



k B T 
J e = D T eV 



D T e\7 



k B T : 
2 1 



— fi e 

T 7n + Ah-V — , 

k B l p k B l 



(56) 
(57) 



where we have omitted the arguments r, t, and have in- 
troduced the particle diffusion coefficient D, the ther- 
mal diffusion coefficient Dt, and the thermal conductiv- 
ity D t h, all of which possess dimensions of length 2 /time. 
The powers of p and e in the prefactors of the gradients in 
(|56|) and (|57j) can be determined from dimensional analy- 
sis, by observing that the left hand side of (l56l) possesses 
units of l/(time x length 2 ) and that of |57|) has units of 
energy/ (time x length 2 ). Note that the "cross terms" , i.e. 
the prefactor of the second gradient in (151))) and of the first 
gradient in ((57)) are identical as requested by the Onsager 
reciprocal relations. See appendix [B] for a derivation of 
(|56|) and (|57p starting from a dissipation function. The 
change in entropy is obtained via the Gibbs-Duhem re- 
lation (l55l) as 



1 



p 



s = — e p 



T 



(58) 
(59) 



where (|59[) follows from the continuity equations (|53|) and 
(El. 

Bearing in mind the structure of the Euler- 
Lagrange equations (|2^|) and (f3"0"|) . we replace p by 
SE[p, s]/Sp(r, t) + v(r, t) and T by SE[p, s]/Ss(r, t). Here 
we have allowed the external potential to be time- 
dependent, in order to model a corresponding external 
influence on the system. Hence we rewrite (|53")) and ([54]) 
as 

D „b a + v D T „ 1 , an , 

(60) 



k B b s 



k B b s 



k B p b s 



(61) 
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where we have used the short-hand notation for the first 
functional derivatives of the internal energy functional, 



SE\p,a] 



5p(r) 
SE\p, s] 



Ss{r) 



p(r,t),s(T,t) 



p(r,i),s(r,t) 



(62) 
(63) 



Performing the replacement of the local temperature and 
the local chemical potential by the corresponding func- 
tional derivatives in (l59l) yields 



(64) 



The equations for the currents and (fBTj) together 
with the continuity equation for the particle density (|53[) 
and for the energy density (f54|) . along with (|64|) . form a 
closed set of equations for the time evolution of p(r, t) 
and s(r,t), for given v(r,t) and initial conditions p(r, 0) 
and s(r, 0) at time t = 0. In general the diffusion coeffi- 
cients D,Dt and D t ^ will depend on p(r,t) and s(r,t); 
assuming them to be constant would be the simplest ap- 
proximation. 



V. CONCLUSIONS 

In summary, we have developed a variational formu- 
lation of classical statistical mechanics, which is cen- 
tered around the internal energy as a functional of the 
one-body density distribution p(r) and the position- 
dependent entropy distribution s(r). Although the defi- 
nition of s(r) is not unique [cf. Eq. (|13p for the equilib- 
rium value so( r )]j the current choice possesses two im- 
portant properties that make it a suitable variable in the 
variational framework: i) the space integral of s(r) is the 
macroscopic entropy, and ii) the definition is local in the 
sense that it probes the entropy under the condition that 
a particle resides at the space point r considered. One of 
the Euler-Lagrange equations for the minimization of the 
grand potential is very similar to that of DFT based on 
the Helmholtz theory, i.e., the functional derivative with 
respect to the density field is related to a local chemical 
potential, cf. (|29p. Physically, such a situation can be 
realized by an external potential acting on the system. 
The internal energy functional E[p, s] depends on the 
local density p(r) and on the entropy distribution s(r). 
The functional derivative with respect to s(r) gives the 
(constant) temperature in equilibrium, cf. (|30[) . Having 
this further Euler-Lagrange equation is to be considered 
a strength of the theory, when it comes to applications 
using an approximate functional. Rather than having to 
implement the physics of T = const on the level of the 
approximation for the free energy functional, the inter- 
nal energy functional offers an additional mechanism to 
relax to equilibrium via an inhomogeneous entropy dis- 
tribution. 



In Levy's constrained search method, which we used 
for formulating the variational framework, there is no 
need for introducing a field that is conjugate to the lo- 
cal entropy distribution. Hence the situation is different 
from the local chemical potential that is conjugate to 
the one-body density. The relationship between these 
thermodynamic fields plays a crucial role in the standard 
Mermin-Evans formulation of DFT. However, in equilib- 
rium there is at least no simple conjugate to the entropy 
distribution. Such a role would be played by an position- 
dependent temperature, which we deliberately avoided in 
the derivation presented in Sec. HU 

Obtaining dynamical equations for the density and en- 
tropy distributions is straightforward when using linear 
irreversible thermodynamics in a continuum description 
as a starting point and replacing the fields for tempera- 
ture and chemical potential by the appropriate functional 
derivatives of the internal energy functional, cf. (|56[) and 
(|57p . This approach is phenomenological and we have 
made no attempts at deriving the dynamics from first 
principles under controlled approximations for the mi- 
croscopic dynamics. While the structure of the dynamic 
equations is a straightforward extension of dynamical 
DFT, there is also an important distinction: When us- 
ing the Helmholtz free energy functional, in principle any 
(non-pathological) density field is a physically realizable 
one via choice of an appropriate external potential. The 
situation is different when considering the internal energy 
functional and prescribing both the density field and the 
entropy field. In general, no corresponding equilibrium 
situation will exist, i.e., one cannot choose an external 
potential and a temperature, cf. (|2"5|) and (|3TJ)) . so that 
the given trial fields p(r) and s(r) become equilibrium 
quantities. This effect is far less subtle than that of rep- 
resentability of trial density fields, cf. the discussion in 
.15]. 

Clearly, true hydrodynamic effects, that originate from 
local momentum and angular momentum conservation, 
are neglected in the treatment of Sec. lIVl However, there 
remains a wide range of interesting physics associated 
entirely with diffusive dynamics in (complex) liquids 
2C|, see e.g. Dhont's treatment of thermodiffusion 
22jj. 



We have formulated a variety of standard DFT approx- 
imations in internal energy language. The mathematical 
structure of some of these functionals appear unfamiliar 
in a variational context, cf. the form of the ideal gas in- 
ternal energy functional and the way in which the ideal 
gas and interaction contributions are coupled in the case 
of hard spheres, cf. (|4"Tj) . Other approximations are con- 
sistent with expectation, i.e., the addition of a mean-field 
energy contribution (145[) . the local density approximation 
(|46l) and the Taylor expansion up to second order around 
a homogeneous (fluid) state (|4"T)) . 

The potential importance of the current work lies i) in 
the additional insights that can be gained from studying 
the entropy distribution in applications within existing 
approximations such as these described in Sec. IIII1 and 
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ii) in the possibility to construct internal energy function- 
al that couple the density and entropy contributions in 
novel ways. Investigating the implications for the dy- 
namical test particle limit 0, [24| is an interesting topic 
for future work, as is considering quenched-annealed mix- 
tures (25l427l | and the dynamics of atomic liquids [H| . Fi- 
nally note that changing the thermodynamical potential 
as we have done here is very different from changing to a 
different ensemble, see e.g. [29| for DFT in the canonical 
ensemble. 

Note added in proof. - The current theory possess simi- 
larities, but also significant differences to the approach by 
Phil Attard [3l|. His theory is, broadly speaking, based 
on the entropy functional with the internal energy being 
a variable. 
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Appendix A: Fundamental-measure theory 

Although Rosenfeld's functional [for a hard sphere 
mixture with one-body density profile Pi{v) of species i] 
possesses the structure of (j40|) . 



where [N(x)] m is the m-th matrix power of 



Fexc[{pi}} = J dx$(x), 



(Al) 



the position coordinate x is very different from the ar- 
gument r of the entropy distribution (|16[) . Rather than 
corresponding to a particle position, x is a mere convo- 
lution integral that couples the FMT weight functions in 
order to represent the hard sphere Mayer bond, and for 
third and higher orders in density, constitutes the cen- 
ter of star diagrams that are formed by weight function 
bonds [H. 

Both the definition of the entropy field (fT6]) . and in 
the ratio of entropy and density distribution in the ex- 
ponential of the hard sphere internal energy functional 
(|41[). one expects a particle to be located at the position 
considered. In FMT there is no shortage of position inte- 
grals over the density, hence the problem is to single out 
one of them in a non-biased, "symmetric" way. 

In order to achieve this we start from the power series 
of FMT HI, which reads 



N(x) 



/n 3 (x) ra 2 (x) rax (x) «-o(x)\ 

n 3 (x) ^ rxx(x) 

n 3 (x) n 2 (x) 

V n 3 (x)/ 



(A3) 



Here the weighted densities n Q (x) are obtained by con- 
volution, 



*a(x) = J dr Pi (r)w^(r - x), 



(A4) 



where the w a (-) are the Kierlik-Rosinberg FMT weight 
functions [7] , and the sum is over all hard sphere species i. 
Note that F exc [{pi}] in (IA2[) is a 4 x 4- matrix and that 
the physically relevant entry -F e xc[{pi}] is that in the first 
row and last column [30j |. 

We rewrite the integrand in (|A2[) as 



N(x)E^ TT^Wr- 1 (A5) 

* — ' mm — 1 

m— 2 v ' 

= £ / drWW (x - r)Mr) £ 7777777"^ N*)]™" 1 , 

i J m=2 V ' 

(A6) 

where the matrix of weight functions is defined as 



^(x) «,W(x) «,<*>(x) W «(x)\ 

«£>(x) ^ 

«4 ( x ) 4 l) ( x ) 

V wf{x)J 



(A7) 



Reintroducing the x-integral and re-arranging in (|A6j) 
gives 

F C xc[{p,:}] = k B T ( dr£ Pl (r) 

^ i 

r 00 1 

x /dxW^x-r)^— -[N(x)]'- 1 , 

(A8) 

= k B T f drJ2pi(r)*i(r), (A9) 

i 

where we have defined the free energy density (per par- 
ticle) for species i as 

*<« = / rfxW 4 (x-r) £ ^-^[N(x)]- 1 . 

(A10) 

Hence we single out the entry in the top right corner of 
the matrix in (IA9I). and rewrite it as 



: [{ Pi }} = k B T fdxf^ J_ [N(xf (A2) F exc [{ Pi }} = I dr£>(r) £ * cf> a ) (r), (All) 
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where the asterisk denotes the convolution, and explicit 
expressions for the </> Q (x) are as follows: 



io = l+f^-lJln(l-n 8 ), 

k n2 n2 W1 ^ 
n 3 n z 3 



(A12) 
(A13) 



[A 



In(l-ns) 



m n|(2-n 3 ) 



no 2riiri2 
772 T3 



47rn| 



rt 3 87rn|(l — rt 3 ) ' 
(A14) 

ln(l - n 3 ) 



no 
n 3 



nin 2 (2-n 3 ) n|(2ri3 - 9n 3 + 6) 
nf(l -n 3 ) 247rdj(l - n 3 ) 2 



(A15) 



In summary, the integrand in (jAllI) forms a suitable 
choice for the desired quantity, i.e., 

SHs([{/>i}],r) = -fc fl *(r) (A16) 

(A17) 

Appendix B: Dissipation function 

As a consistency check on (f5tj)) and (|57|) [and hence 
([60]) and (jSTj) ], we derive the currents from a (scaled) 
dissipation function R [l?], [H| , which we assume to be 
given by 



2p 



(Bl) 



Here i? is a scaled object with dimensions of (length x 
time) -1 . One can verify explicitly that (|B1|) generates 
the expressions ([56)) and ([57]) via 



OR 



OR 



d 



(B2) 
(B3) 



Furthermore, one can show explicitly that for the entropy 
production 



s + V • J, = 2k B R, 



holds, where the entropy current is 



(B4) 



T 



(B5) 



consistent with the Gibbs-Duhem relation 



Finally note that and (|57|) can be written as a 
matrix product 



K\ (Dp D T e \ /V^ 



(B6) 



where the matrix of kinetic cocfhcicnts [on the right hand 
side of (|B6[) ] is symmetric, as requested by the Onsager 
reciprocal relations [13, [3 ■ 
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